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Collimated relativistic outflows, or jets, are amongst the most energetic 
and relativistic phenomena in the Universe. They are associated with su- 
permassive black holes in distant active galactic nuclei (AGN) [lj, accreting 
black holes and neutron stars in binary systems [2j and are believed to be 
responsible for gamma-ray bursts |3j. The physics of these jets, however, 
remains something of a mystery in that their bulk velocities, compositions 
and energetics remain poorly determined. Here we report the discovery of 
an ultra-relativistic outflow from a binary accreting neutron star accreting 
gas within a binary system. The velocity of the outflow is comparable to the 
fastest-moving flows observed from active galactic nuclei, and its strength 
is modulated by the rate of accretion of material onto the neutron star. 
Shocks are energized further downstream in the flow, which are themselves 
moving with mildly relativistic bulk velocities and are the sites of the ob- 
served synchrotron emission from the jet. We conclude that the generation 
of such highly relativistic outflows does not require properties unique to 
black holes, such as an event horizon. 



Circinus X-l is a bright and highly variable X-ray source, containing a stellar-mass 
compact object accreting from a binary companion star. X-ray dips and outbursts 
with a period of 16.6 days are most readily interpreted as enhanced accretion during 
periastron passage of the accreting object in a significantly elliptical binary orbit |H|3]. 
Observations of type-I X-ray bursts indicate that the accreting object is a neutron 
star 0. In the 1970s Cir X-l was also associated with bright radio outbursts [7JIH]. 
Since then, there has been a decline in the strength of the radio counterpart, which has 
been found to reside within an extended radio nebula 0. Structures on arcminute- 
scales have been imaged within this nebula jTUj, and more recently a one-sided jet on 
arcsecond-scales, which aligns with the larger structures, has been discovered [TTj . 

Bright radio flares associated with the production of a relativistic jet are now es- 
tablished to be a common property of accreting black holes in outburst, but are less 
commonly associated with accreting neutron stars jS] . In several cases the radio emis- 
sion has been resolved into outflowing components, sometimes either side of a stationary 
core, with bulk Lorentz factors (the Lorentz factor T = (1 — /3 2 ) -1 / 2 , where velocity 
v = (3c and c is the speed of light) typically inferred to be in the range 2 < V < 5 
[2|. Based on these observations a generic qualitative picture has emerged in which 
electrons are accelerated to relativistic energies close to the accreting compact object, 
and then ejected from the system in a bipolar flow, from which we observe synchrotron 



emission until energy losses, primarily due to expansion of the plasma, reduce the 
emission below observable levels. However, very long baseline interferometric obser- 
vations of the neutron-star binary Sco X-l [T2*l IT3] have revealed knots with mildly 
relativistic bulk velocities {(3 ~ 0.5) which are energised by a beam which is itself not 
directly detectable, but must have a Lorentz factor T > 3. Furthermore, it has been 
suggested that the observed synchrotron emission results from a moving shock (as 
discussed for AGN and Gamma-ray bursts), a possibility supported by recent observa- 
tions of large-scale X-ray jets from a transient black hole X-ray binary, revealing in-situ 
acceleration of electrons to TeV energies [To] . 

At all frequencies and angular scales observed, Cir X-l is associated with an ex- 
tended radio structure. Careful inspection of Figs 1 and 2 reveals that X-ray flares are 
associated with a brightening of the 'core' of the extended radio structure. Previous 
comparison of Hubble Space Telescope and radio images ^JJ confirms the association 
of the radio 'core' with the optical counterpart of the binary, and so we associate the 
extended emission/knots with an outflow from the system. The radio spectral indices 
(—1.1 < a < —0.7) indicate optically thin synchrotron emission. 

We shall focus initially on the observations of October 2000 (left-most panel of Figs 
1 and 2), which have a sampling of every second day. In the 'quiescent' observation 
prior to the X-ray flaring the radio source at 4.8 GHz is clearly associated with a 



bright knot near to the core (best radio position from plus a weaker extended 

structure. At 8.6 GHz the extended emission is resolved into two components, knots 
A and B, which are not co-aligned with the core. Following the X-ray flaring, the 
core has brightened significantly within one day, and within three days the extended 
emission is brightening, peaking within five days of the outburst. Subsequently the 
extended emission fades. In May 2001 (middle panel of Fig 1) we performed a similar 
set of observations with two-day sampling and saw a similar effect. Following the X- 
ray flaring the extended emission brightens at 4.8 GHz within two or three days, and 
there seems to be a 'bridge' of radio emission from the core to the knot immediately 
after the flare; in addition in the difference maps (Fig 2) there is some indication of 
a counterjet. In Dec 2002 we made a sequence of observations with daily sampling. 
This time following the X-ray flare the core brightened strongly and a resolved bridge 
of radio emission appeared along the apparently curved axis of the extended emission. 

At all epochs we clearly see evidence for the transmission of energy from the binary 
core to the extended radio emission following an outburst. To what apparent velocity 
does this signal correspond ? For a proper motion /z, the apparent velocity expressed 
as a fraction of the speed of light /9 app = t>a PP /c, for the source at a distance d is 



All the maps show the brightening of the extended components within days of the 




X-ray flare. At all epochs the extended emission at 4.8 GHz, at a distance of 2.2-2.5 

arcsec from the core, brightens within three to five days. From these numbers we can 

estimate a conservative lower limit to the proper motion associated with the energising 

beam of 400 mas d _1 . The most likely distance to Cir X-l is ~ 6.5 kpc ^U]. From this 

we can establish that the apparent signal propagation velocity between the core and 

knot is /9 app > 15. For an intrinsic velocity /3 at an angle 9 to the line of sight, 

/3sinfl 
^ app ~ (l-/?cos0) 

which has a maximum value of /3 app = T(3 when f3 = cos 9. Since the observations 
indicate /3 app > 15, this implies Yj3 > 15, for which the minimum velocity solution is 
r > 15 (i.e. f3 > 0.998). This is by far the most relativistic flow observed to date within 
our Galaxy. Furthermore, for any value of (3, solutions may only be obtained for angles 
to the line of sight 9 < 5°, indicating that the signal responsible for the brightening 
of the knots must an angle very close to the line of sight. Furthermore, 

the images, in particular at 8.6 GHz, give a clear indication of bending of the jet; by 
comparison with AGN we can also infer a small angle to the line of sight, in which case 
projection effects increase the apparent bend angle. 

This lower limit tothe Lorentz factor may correspond to the bulk velocity of the 
flow between the core and the jets, or of a pattern (shock) propagating along this flow. 
In the latter case, in the approximation of a ID relativistic shock propagating along 



a relativistic fluid, the bulk (Tbuik) and shock (r s h oc k) Lorentz factors are related by 
Tbuik = r s hock/v^2 in the limit (appropriate here) of /3 s h oc k - > 1 (e.g. [TOJ ) . Thus our 
lower limit to the signal Lorentz factor of T > 15 may place a lower limit on the bulk 
Lorentz factor of the flow of Tbuik > 10. 

Further clues to the nature of the sites of emission are provided by linear polarisa- 
tion maps, to be presented elsewhere, which reveal a magnetic field in the knots which 
is parallel to the jet axis in the core and perpendicular to the jet axis at the knots. 
This supports an interpretation of the extended emission as synchrotron emission from 
shocks formed at a longitudinal compression of the jet flow, as in AGN ^7]. Further- 
more, it is clear that there is secular evolution of the mean radio structure between the 
epochs, with a general trend for the propagation of components away from the binary 
core. 

These observations reveal that underlying the observed radio jets of the accreting 
neutron star Cir X-l is a flow with a bulk Lorentz factor > 10. In another neutron 
star binary, Sco X-l the underlying physics may be similar |T21 ITS] . Sco X-l is the 
prototype of the class of 'Z sources', representing the six most luminous neutron star 
accretors in our Galaxy, and it has further been suggested that Cir X-l itself shares 
some of the properties of this class [TH]. Thus it seems plausible that production of 
such ultrarelativistic jets may be a generic feature of accretion onto neutron stars at 



or near to the Eddington limit - this is at odds with the popular suggestion that the 
velocities of jets from neutron stars (and indeed all accreting objects) would be limited 
to the surface escape speed, in this case ~ 0.3c [T3] . 

The minimum energy associated with brightening of the extended emission is E m m ~ 
10 40 erg for a source volume of radius one light day (the estimated energisation timescale); 
the associated power input rate is Pj > 10 35 erg s -1 . If the unseen energising signal 
were isotropic - for example a spherical blast wave from the binary - this would cor- 
respond to a total power output rate of > 7 x 10 39 erg s -1 , more than an order of 
magnitude above the Eddington limit for a neutron star. Based on these arguments 
we interpret the signal as arising in a collimated ultrarelativistic flow. 

The velocities of jets from black holes in our galaxy are not well determined; in 
particular it is almost impossible to place an upper limit on the Lorentz factors of the 
observed motions [211] ■ In the case of AGN it has been argued that the bulk Lorentz 
factor of motion is T ~ 10 , although in a small fraction of blazars apparent velocities 
(3 > 30/i -1 (where the Hubble constant H = lOO/i km s _1 Mpc _1 ) have been measured 
on VLBI scales with the same implications for the underlying bulk velocity as 
derived here for Cir X-l. Thus these observations demonstrate that the study of jets 
from galactic XRBs, whether containing neutron star or black hole accretors, may be 
directly applied to our study of jets from the most powerful engines in the universe, 



the AGN. Furthermore, the discovery of such ultrarelativistic flows associated with a 
neutron star points us clearly to the common features of neutron star and black hole 
accretion, e.g. the accretion flow, and not properties unique to black holes, e.g. an 
event horizon or ergosphere, as the sites of the key physics in jet formation. 
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Fig 1 An ultrarelativistic outflow: sequences of radio observations of Circinus X-l in October 
2000, May 2001 and December 2002. At each epoch, observations were made simultaneously 
at 4.8 and 8.6 GHz. White tickmarks indicate time steps of one day; the blue bar indicates the 
time of the X-ray flaring as observed by the Rossi X-ray Timing Explorer All-Sky Monitor. In 
October 2000 and May 2001 the observations were spaced every two days; in December 2002 
they are daily. At each epoch the u-v coverage of the radio observations is identical for each 
image; maps in October 2000 and May 2001 are 'uniformly weighted', those in December 
2002 are 'naturally weighted'. The crosses indicate the location of the binary 'core' from 
jllj . and their size is proportional to the 'core' radio flux density. The observations reveal 
that following the X-ray flaring the extended radio structure brightens on timescales of days. 
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Fig 2 The same data as Fig 1, but presented in the form of 'difference' maps. At each of the 
three epochs the top (reference) image at each frequency is the same as that presented in Fig 
1. The contour images below it correspond to the maps from the corresponding frames in 
Fig 1 with the reference image subtracted, in order to clearly illustrate changes in the radio 
structure. Solid lines indicate positive residuals, dashed lines indicate negative residuals; 
contour levels have been chosen to indicate the level of noise in each image. At each epoch, 
the strongest variability in the difference maps is clearly occurring in the region of the jet 
/ knots. Solid ellipses in the lower right corners of the bottom panels indicate the size and 
shape of the synthesised beam at each epoch and frequency. 



